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Applying the bosonization procedure to weakly coupled Hubbard chains we 

discuss the fixed points of the renormalization group procedure where all spin 

G , excitations are gapful and a singlet pairing becomes the dominant instability. 
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I. INTRODUCTION 

Recently the field of non-Landau Fermi liquid states in various quasi-one dimensional 
systems has been very active. Although basic properties of purely one dimensional systems 
(chains) are quite well known by now, it still remains to be understood how these change 
under coupling between chains. In particular, if an infinite number of chains is coupled to 
form a two dimensional array then some kind of a dimensional crossover occurs. Despite 
numerous intensive studies of these questions, it still remains open how the non-Landau 
Fermi liquid one dimensional features evolve to isotropic two dimensional behavior. 

On the other hand one may expect new physics in a system with a finite number of cou- 
pled chains which may exhibit an unusual amalgamation of both one- and two-dimensional 
features. Besides a general theoretical interest, such systems are also attractive because 
they can be found in real materials. Recently it was pointed out that some substances such 
as Sr2CuiOQ provide a physical realization of weakly coupled double chains [|l|. Moreover 
higher stoichiometric compounds in the series Srn-iCun+i02n present examples of coupled 
iV- chain ladders. 

This expectation is reinforced by the behavior of S* = 1/2 Heisenberg multichain or ladder 
systems. Whereas the single spin chain shows quasi-long range order with gapless spinon 
excitations, the double chain system shows spin liquid behavior [^, |^, ^ with strictly 
short range order and finite gap in the spin excitation spectrum. This contrast in behavior 
has led to the conjecture that a lightly doped double-chain system should preserve the spin 
gap and become superconducting [jl| . Further in view of the insensitivity of the spin liquid 
state to the ratio of interchain to intrachain coupling [Q, one expects that it is robust and 
should occur also for a Hubbard model even in the weak coupling regime. This expectation 
is strongly supported by a recent numerical study [Q. This is our motivation to examine the 
renormalization group (RG) theory of weakly coupled Hubbard double chains and to look 
for a spin liquid fixed point. 

Recent weak couphng RG studies of the double-chain Hubbard model H revealed some 



strong coupling fixed points characterized by enhanced singlet pairing. However the anal- 
ysis performed in ||^ is essentially restricted by the case of weak interactions between the 
fermions. These authors did not examine the half-filled case with weak interchain hopping 
where the Umklapp processes on the individual chains become relevant. 

In the present paper we undertake an attempt to construct a description of the spin gap 
fixed point by using a bosonic representation. This is expected to be an adequate tool to 
demonstrate a development of the strong coupling regime in both cases of weakly coupled 
Hubbard chains and of a strongly correlated double chain t — J model. The latter case will 
be considered in a later publication 0. 

To clarify the essence of the double chain physics in the presence of strong correlations 
it is worthwhile to begin with a review of well known properties of the single chain Hubbard 
model. 

Away from half filling the model can be only found in the so-called Tomonaga- Luttinger 
(TL) regime which corresponds to both gapless spin and charge excitations |p]. It is cus- 
tomary to describe the TL behavior in terms of spin and charge correlation exponents Kg 
and Kc- 

The spin exponent Kg equals to unity everywhere in the TL regime while K^. gradually 
increases from the value K^ = 1/2 at values of the onsite repulsion U = oo and any electron 
density p 7^ 1 (as well as at p ^ 1 and arbitrary U/t (t: intrachain hopping)) as U increases 
or p gets smaller. 

In the regime of strong correlation at p close to unity one can argue that gapless spin 
fluctuations drive the coupling constants of the charge sector to the repulsive region {Kc < 1). 
In this case we can expect a change of behavior when the exchange coupling between chains 
is introduced. The spin gap, which we have argued is a robust feature of the single rung 
ladder at p = 1, acts to cutoff the spin fluctuation spectrum at low energies so that these 
may not renormalize the charge couplings significantly. As a result there can be an effective 
attraction at p < 1 without a threshold value of the ratio U/t. This will manifest itself 
in a finite spin gap also at p < 1 and a scaling to the Luther-Emery line rather than the 



Tomonaga-Luttinger line. Such behavior has been reported for single chains models with 
a spin gap caused by frustrated or modulated exchange couplings on the single chain 0, 
p!0| . In the present case of a single rung ladder or double chain, we will choose the regime 
with the interchain hopping t_L weak {t± << U) so that as p ^ 1 real interchain kinetic 
energy processes characterized by tj_ will scale to zero. In this limit the induced interchain 
exchange processes will remain finite and dominate. 

A mean field analysis |jll| of the t — J model predicts that the spin gap remains upon 



doping and the RVB state at p = 1 evolves into a superconductor with approximate d-wave 
symmetry. Noack et al [^ found a similar behavior in their numerical studies of moderately 
coupled Hubbard ladders. 

In this paper we investigate the case of a weakly coupled Hubbard ladder with p ^ 1 
and t± < U < t using RG methods and look for a strong coupling fixed point with the same 
characteristics. 

II. BOSONIZED WEAK COUPLING LIMIT OF THE DOUBLE-CHAIN 

PROBLEM 

To get the first insight into the problem we start with a conventional bosonization of the 
small U/t Hubbard model on two weakly coupled chains. 

+UY^{u\^Ui^u{_^u,^^^ + d\ji^d{_Ji^.^) (2.1) 

Here Uia and di^j denote fermions on upper ("u") and lower ("d") chains. 

Apparently, at f/ << t_|_, t the interaction term has to be treated as a small perturbation 
to the rest of the Hamiltonian (2.1) and the bare transverse hopping leads to the formation 
of two (bonding (B) and antibonding (A)) bands: {A^B) = 4=(u =p (i). Thus at U « t± 
a proper starting point is provided by the two-band model which was previously studied in 
the framework of a general weak coupling "g-ology" |T^, ||T3[. The analysis carried out in 
M was also based on the two-band picture. 



However, if the opposite condition t± << U is satisfied then the effect of band sphtting 
is completely suppressed due to the requirement to avoid a double on-site occupancy. This 
behavior persists down to the quarter filling ( p = 1/2). In the framework of the RG 
approach this phenomenon manifests itself as a vanishing of the renormalized t±. In view of 
this we suppose that in the case t± « U one has to start from the picture of two degenerate 
bands to implement correctly the fermion correlations. 

The preceding RG analysis of the general two-band model in absence of Umklapp pro- 
cesses |12|, |T^ already shows many technical complexities. For this reason the results of 



these studies are not simply physically transparent. Moreover it turns out that all nontrivial 
fixed points are located far in the strong coupling regime where the lowest orders RG calcu- 
lations cease to be valid. So one might expect that a more informative investigation can be 
done on the basis of a bosonic representation which is usually capable of yield of a correct 
evolution toward strong coupling and even to provide an asymptotically exact solution of 
the Luther- Emery type ^^. Recently the method of bosonization was applied to the double 



chain problem in the context of a special model which includes only forward scattering []T3 . 



This analysis |T5| led to the prediction that coupled chains provide a proper basis for the 



occurence of singlet pairing. 



In this paper we will perform a more general analysis than that of [^ to see whether 
the above statement holds for a wider class of models. 

To proceed with a bosonic representation we introduce a conventional set of bosonic 
fields 0{, 0{ where the "flavor index" / has one of two values u or d. These fields describe 
fluctuations of charge (c) and spin (s) densities respectively. In the continuum limit the 
fermion operators can be written in terms of these variables as follows 



z 



^i{x) ~ ^ expizi^kpx + ^(/i0{ + ef + fia(f){ + aO^)) (2.2) 

where 9^^ are the dual fields to the (j){.g {d^6{g = e^udv(i){s). 

Applying the formula (2.2) and introducing the hnear combinations 0^^ = -75 (0"* ^ 
Of,. = 75(6'",, ± 0^,.) corresponding to total (" + ") and relative (" — ") charge or spin 
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6^ , , 



density fluctuations, one can readily obtain the bosonic form of the Hamiltonian (2.1) 

+t_L(cos0^ cos0~ + cos(0^ + 5x) COS0J") cos6'^ cos6'g" + 

qbs cos 20^ cos 20J + gu cos(20^ + 25x) cos 20~ (2.3) 

where last two cosine term represents spin backscattering and Umklapp processes respec- 
tively. Each of these terms is, in fact, a sum of two contributions cos2v20cs coming from 
u and d species. The Umklapp term becomes relevant when the doping 5{= ^ — kp) van- 
ishes. As usual, the bare values of the correlation exponents Kf^ = ^ can be changed 
by shortwavelength renormalizations. Neglecting these corrections we obtain that the bare 
correlation exponents governing the charge dynamics K^{= (1 + ^)~^^'^) are smaller than 
unity while the spin exponents Kf{= (1 + ^Y^"^) are opposite and Kf > 1. In addition, 
the bare values of the amplitudes gss and gu are equal to ^. 

To perform a renormalization procedure we divide up all variables on slow and fast 
components and then integrate out the fast variables. Using the bare values of correlation 
exponents one can estimate scaling dimensions of various terms in (2.3) according to the 
conventional formula [|TB|(7 and 7' are arbitrary): 

A(cos70t cosYC) = ji^'^t + t±) (2-4) 

c,s 

The Hamiltonian (2.3) has to be supplemented by extra terms which are generated in the 
course of renormalization. Indeed, performing an expansion of the partition function Z = 
Trexpl—PHs) in t± one immediately observes that the interchain hopping produces the 
following relevant terms (new couplings gi should not be confused with the traditional g- 
ological notations): 

AH = gi cos 20~ cos 26~ + (72 cos 26~ cos 20~ + g^ cos 26~ cos 20j" + 
g^ cos (20^ + 25x) cos 26^7 + g^ cos 20^ cos 20 J + g^ cos 20^ cos 26*7 + 
g-i cos(20+ + 25x) cos 20' + 5-8 cos(20+ + 25x) cos 20+ + 5(9 cos 29' cos 29' (2.5) 



All these terms have scaling dimensions not greater than two and result from the second order 
perturbation corrections to the single chain Hamiltonian AH ~ t^ < {Y,i Uidi + dlui)"^ >. 

Physically these terms correspond to processes of coherent interchain particle-hole and 
particle-particle hopping triggered by the single particle one. The crucial importance of 
these processes was previously pointed out by many authors (see, for instance, |T^, [|T^]). 

In the second order in t± the RG equations derived by the use of the method of [|l9l have 
the following form (^ = Inx): 

I = (2 - /C - ^)«, + |(Ar + ^ - K; - ±-_) - ,,«„ (2.6) 

I = (2 - /C - ^),, + |(/C + ^ - /C - ^) - ,,,ss (2.7) 

f = (2 - K - ^),3 + |(/C + ^ - K- - ^) - J.9BS (2.8) 

f ^ (2 - Kt - ^)«. + |(/C + ^ - Kt - ^) - ',.,. (2.9) 

t - p - '^'= - '^'- )«= 4'*^ + *• - 7^ - 1?F' <""* 

^ = (2 - A-+ - -i^)j„ + |(A7 + -i^ - A-<+ - ;^) - Ms - ftS, (2.11) 

^ ^ (2 - A- - ^)«. + |(A- + ^ - /C - ^) - 9.* (2,12) 



,^ = (2 - /C - ^)9> + |(/C - ^ + /C - 7^) - m,. (2.13) 



(2.14) 



^ = (2 - ir; - Kt)gBs - 929, (2.15) 

{2 - K- - Kt)9u - 9i9^ (2.16) 



d9u 
di 



^^^ = \{-K-{9l + 9l + 9l) + ^{9l + 9l + 9'r + 9l)) (2.17) 

^^^ = \{-K;{9l + 9l + 9ls) + ^{9l + 9l + 9l + 9l)) (2.18) 

-Kt{9l + 9^7 + 9l + 9u) (2.19) 



d\0gK+ 1 ^^+2 ,2,2,2- 



di 2 

^^^ = -\Kt{9l + 9l + 9l + 9ls) (2.20) 

^^2-i(.,.^..;.^) (2.21) 

In addition, there are two equations describing evolutions of velocities Vc^g but one can always 
include these corrections into the definition of the correlation exponents. 

In comparison with the equations obtained in |T2[, |T^ our RG equations (2.6-2.21) 



are already written in terms of physically relevant combinations of original "g-ological" 
couplings, so one could hope that this description might appear to be more transparent. 
As we shall show the above system consistently demonstrates a development of the strong 
coupling regime in rather general conditions, so we don't account next-to-leading order 
corrections which can be only necessary if one discusses fixed points at finite coupling. 

First, in the case of spinless fermions away from half filling the only relevant couplings 
are 91,92 and K~ and the equations (2.6-2.21) reduce to those found previously pO], pT|, 



22 



However in the physically relevant case of spin one half fermions away from p = 1 one 
can only neglect couplings g^, gj, gg and gu associated with Umklapp processes and then the 
number of residual couplings is large (10) and coincides with that found in ^. 
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The fact that eq. (2.6-2.21) originate from the repulsive Hubbard model simplifies their 
analysis significantly. To see that one can choose a two-step renormalization procedure to 
that of Ref. |^ and integrate the above equations first up to the scale ^o = In ^ where 
the renormalized amplitude of single particle hopping t±{^) becomes of order unity (and 
stops). It can be easily seen that at ,^ ~ ,^0 one can still neglect renormalizations of the 
correlation exponents K^ = 1 — ^ and Kf = 1 — ^ from their bare values corresponding 



to zf = -zf = A = -^. 



By straightforward generalization of the analysis of |^ , ||2^ one can obtain the evolution 



of couplings gi{C) given by the eq.(2.6-2.14) with only inhomogeneous terms proportional to 
t]_ kept 

2At^5 _ A,? 

^.(O = ati(0)^^^--^ (2.22) 

where Aj, i = 1, ..., 9, BS, U denote dimensions of relevant operators and Ci are the coeffi- 
cients standing in front of terms proportional to t]_ in the r.h.s. of (2.6-2.14). 

It follows from (2.22) that (71(^0) = ~5'2(^o) = ~fi'3(^o) = fi'4(^o) = —A/2 while all the 
other couplings gi{^Q)-,i = 5,..., 9 are of order A^. On the other hand at ^ > ^0 one can 
also omit in (2.6-2.14) all inhomogeneous terms using g{^o) as bare values. Naively, it could 
mean that one has to account the leading couplings fifi, 2,3,4 plus gss^Qu first and then to 
treat all the rest as additional perturbations. However it turns out that the solution is not 
so straightforward. 

Let us consider first the case away from half filling. Then it can be shown that couplings 
g2,gz and gss all tend to zero though gi diverges. Asymptotically the following relations 
hold 

^S^ ~ exp({(A-- + -L - K; - K:)) ^ 
9iK) K. 



' - Ki()- ''B:'''°\ (2.23) 



Thus we infer that K~{^) vanishes while K~{C,) goes to infinity. But this means that the 
assumption about smallness of couphngs (75,. ..,9 made on the basis of their values at ,^ = ,^0 
was not quite correct. Namely one has to include those terms which contain one of the 
fields 0^ or 6 J which are "close" to get locked. A simple inspection yields that the second 
relevant coupling (besides gi) is Qq while in the case of p = 1 one has to keep g4^8,u ^^ well. 
The resulting system of equations in the range C,o < C, < j reads as 

-^ = -^i^c - ^7)91 - 9i9u (2.24) 



(2.25) 



dge 
d^ 


4(. 


. — Z^ 


,)9&- 


-9i98 




dz~ 
d^ 


= 9\ 


+ 9u 




1 


dz- 
di 


-9\- 


-9l- 


9l 




dzt 
di 


= 9l 


+ 9I 






dzt 

AC 


9l + 


9I + : 


9l 



(2.26) 



(2.27) 



f2.28) 



(2.29) 

Away from half filling all (74,8,(7 freeze at .^ ~ In | (and, consequently, z^ is frozen too) and 
there are only qi^q left over. Then the system (2.24-2.29) demonstrates a development of the 
strong coupling regime in all channels except the "+" charge one (namely, gi{C) ■, 9&{C) ~^ ~oo 
and -2(7(0) -^^(0 ~^ 00 while z~{^) -^ —00). As usual, these tendencies have to be undersood 
in such a way that at ,^ ~ |- all couplings reach values of order unity and don't vary further. 
Including couplings g^^^u at 5 ^ one can see that they don't alter the behavior found 
for the doped case while the "+" charge sector is also driven to the strong coupling regime 
in accordance with the complete freezing of charge degrees of freedom at p = 1. 
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To facilitate the analysis of leading instabilities of the complete Hamiltonian (2.3,2.5) 
one has to consider eight relevant order parameters where plus and minus correspond to 
intra- versus inter-chain type of ordering 

CDW+ = Y: ^it^^^,. ~ cos(0+ + re) cos(0+ + 0;) (2.30) 

CDw. = y: ^il^zi^ ~ cos(0+ + e-) cos(0+ + e;) (2.31) 

SDW+ = Y ^il^^M,- ~ co<</'c^ + cos(e + e:) (2.32) 

SDW^ = Y ^it^ll-. ~ cos(0+ + 9-) cos(e + 0;) (2.33) 

SS^ = Y: cT^ia^L,,-. ~ cos(e + 0;) sin(0+ + 0;) (2.34) 

SS^ = Y^H^'^-i-<^ ~ cos(e + 0-) M<l>t + 07) (2.35) 

TS+ = Y ^^i.^^M,<x ~ cos(e + ^;) sin(e + 9;) (2.36) 

T5_ = Y ^^i.^-/.,. ~ cos(e + 0c") sin(e + 0;) (2.37) 

We remind that if any of the fields (0^^ or 9^^) gets locked then the corresponding cosine 
acquires a nonzero expectation value and < cos 0(a;) cos 0(0) >^ | < cos 0(0) > ^ as a; 
tends to infinity. On the other hand, fluctuations of both this variable and its dual one 
become gapful. Formally one can identify the state where 0^^ is ordered with the limit 
K^g -^ while 9^^ becomes ordered at K^^ -^ oo. 

Then one can easily see that in the case when 9j and 0^ are locked the only competing 
instabilities are the interchain charge density wave corresponding to the order parameter 
CDW- ~ cos(0^ + 9'^) and the interchain singlet pairing described by SS- ~ cos(^^ + ^jT)- 
In fact, the former state can be also recognized as a counterpart of the two-dimensional flux 
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phase. Indeed, this state is characterised by the commensurate with density flux $ = 2kF 
which is defined as a circulation of a phase of the on-rung order parameter < uldi + rfjuj > 
through a plaquette formed by two adjacent rungs of the ladder. In the case of spinless 
fermions this type of ordering called " Orbital Antiferromagnet" was first considered in ||2^ 
as a prototype of recently proposed two-dimensional fiux states. 

Although the fiux phase can be in principal realized in some extended models we see 
that in our case of the double chain Hubbard model where the field 0~ also gets locked the 
ground state is a spin gapped singlet superconductor. 

It is also instructive to express the above order parameters in terms of the hybridized one- 
particle states corresponding to the abovementioned "bonding" and " antibonding" bands 

CDW. = Y.Ai^AL^-Bl,^BL^ 

cr 

SS. = Y: ^L^l-. - <4,-. (2.38) 

a 

Considering the distribution of signs of the order parameter SS^ on the "four-point Fermi 
surface" {k = {kp, 0), {—kp, 0), {kp, vr), {—kp, vr)) we observe that it corresponds to the "d- 
wave" type pairing. We conjecture that in a two-dimensional array of weakly coupled double 
chains with a continuum Fermi surface this type of ordering does transform into an ordinary 
d-wave pairing. 

III. CONCLUSIONS 

In the present paper we applied the bosonization method to find further arguments in 
support of the recently proposed scenario of singlet superconductivity in the spin gap state 
of doubled Luttinger chains. Previous results obtained in the framework of the mean field 
approach [^ as well as earlier numerical studies [^ , also testify in favor of this picture. 

We also want to stress that our conclusions contradict with a recent claim about the 
existence of the strong coupling fixed point where some spin excitations remain gapless 
made in [03. These authors considered the double chain t — J model without an interchain 
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spin exchange {J± = 0). Then on the bare level their Hamiltonian can be assigned to the 
universality class of the purely forward scattering model considered in |T^ . In this special 
case indeed the only field becoming massive is 6j . In principal, it can't be ruled out that 
for some specific double chain models only part of all relevant fields acquire masses and 
the others remain massless. One example of this type was discussed by the authors of p3[ 



who found only 0^ and (pf to be massive in the framework of the model including solely an 
interchain interaction of fermions with opposite spins. 

However our investigation of the Hubbard-type models shows that the presence of the 
interchain one particle hopping is already sufficient to generate the antiferromagnetic spin 
exchange term with Jj_ ~ .rnaxft u) "^^ich makes all spin modes gapful at p < 1. We believe 
that spin liquid behavior with a finite spin gap which evolves into "d-wave" pairing upon 
doping is a robust feature of a whole variety of isotropic spin models of weakly coupled and 
strongly correlated fermions on double chains. 
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